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■ INTRODUCTION 

\ 

The  statistical  properties  of  signals  reflected  fron  the  ocean 
surface  depends  on  the  stochastic  properties  of  the  surface. 

An  example  of  this  is  the  prediction  of  the  average  energy  reflected 
from  the  surface  as  a function  of  angle  of  incidence.  To  evaluate  this 
quantity,  it  is  necessary  to  know  the  spatial  correlation  of 
surface  heights.  Measurements  or  assumptions  of  the  directional 
wave  spectrum  of  the  ocean  are  first  obtained.  The  spatial- 
temporal  correlation  of  the  ocean  surface  is  expressed  in  terms 
of  this  directional  wave  spectrum.  This  memorandum  describes  a 
method  of  evaluating  this  expression  by  means  of  the  Fast  Fourier 
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'Transform  (FFT).  Since  the  FFT  method  is  much  faster  than  the  usual 
numerical  integration  method,  it  is  highly  desirable  for  experimental 
data.  This  memorandum  presents  the  first  known  evaluation  of 
surface  spatial  correlation  by  the  FFT  method.  _ 

DIRECTIONAL  WAVE  SPECTRUM 


A review  of  the  statistical  properties  of  the  ocean  surface 
is  given  in  References  (l),  (2)  and  (3)  • The  review  will  not  be 
repeated  here.  The  directional  wave  spectrum  is  expressed 
as  A A^C^j  is  the  amount  of  surface  wave  energy 

of  angular  frequency  to  flowing  in  the©  direction  in  the  interval 
d u>d&  • The  relation  between  the  spatial-temporal  correlation 
and  the  directional  wave  spectrum  is 

rr  °° 

r — t Jde>Jdco  e)coAf  ^cose  v-  e ~ co'yj 

^ - n o v d 

where  ^ is  the  acceleration  of  gravity  constant 

u is  the  separation  of  the  surface  points  in  the 
x direction 

V is  the  separation  of  the  surface  points  in  the 
y direction 

'T'  is  a separation  of  tine  at  the  two  points 
(?  is  the  correlation 

S-*-  is  the  mean  square  height  of  the  surface 


Expressing  the  directional  wave  -spectrum  in  terms  of  the  surface 
wave  numbers  k*  feu  in  the  x and  y directions,  we  have  (see  References 
(1),  (2)  and  (3))  d 


€(u)  vj '*')  - k*  <1  k^  k^c-os [k* u +- Vj 


The  spatial  correlation  of  interest  in  this  study  is  obtained 
by  setting  7^0.  Then 

°o  ^ 

<P(Oj  V,  o)  = — ffj  cei  E**^ 

In  order  to  put  this. in  the  more  convenient  Fourier  Transform,  we  have 

*nv‘j  = -L_  laf/Sd  t.  j tg  ax 

^exrC‘t*tTkx»'+  *trk,v'j}  _*y 
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The  expression^  within  the  brackets  is  the  two  dimensional  Fourier 
transform  of  stands  f°r  the  rea^  Pai*t  of. 

Let  us  now  consider  a wind  in  the  x direction.  It  is  sometimes 
assumed  that  the  directional  wave  spectrum  exists  in  the  first 
and  fourth  quadrants  only  and  that  the  spectrum  is  symmetric  about 
the  x axis.  The  above  assumptions  are  intuitively  reasonable. 


<*> 


For  these  assumptions,  ^ 

ei*nv\*rr\/)o')  - i Re  I J1  fd^x  A't** 


")cxf°]Lt  * n 1 + v 


or  r a )x-7 

C(a.  n u1 t 2. IT  o ) - 3-  fte  £ fc*  A’Yfc*,  x n -t-  j 

+ / •»*}  /7t«  A*C  t.,  ^ e*r»F< 1,r  c v'l!  ] 

Let  us  now  consider  the  second  double  integral  and  let 

Thea 

fit*  k/'  + 

- J”  /'c ^7)<SKfoj^C7.TT  C fcy,V  )J 

where  we  have  used  the  symmetry  property  of  A 


Let  F(  ° 
Then 

CC^tTV'j 


- J* /civ;*  AxCk*,  ky)exr°  JjIi.t  (txi/ 

o)=  Re + 


* £~ 


Note  that  Pf',1wjloV)ts  symmetric  in  W'1  but  not  in  general  symmetric 
in  y'  . 

For  many  directional  wave  spectra,  such  as  the  Neumann-Pierson 
spectrum,  Equation  (5)  can  not  be  integrated  analytically.  It 
can  by  evaluated  by  digital  methods f 
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DIGITAL  FORM 

L.et  us  first  evaluate  ft An  digital  form.  The  properties 
of  the  discrete  Fourier  transform  (DFT)  answer  will  contain 
fC . We  choose  a distance  D large  enough  so  that  to  a good 
approximation, 

an  9 

Since  Equation  (6)  is  now  defined  in  a finite  interval,  it  may- 
be approximated  by  a double  finite  stun.  Consider  the  grid 
shown  in  Figure  1. 


o > ^ 
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FIGURE  I 


We  approximate  Equation  (6)  by  a sura  of  rectangular  complex  volumes. 
The  grid  size  h is  chosen  such  that  there  is  not  a significant 
change  of  the  complex  height  in  this  interval.  The  number  of  grid 
intervals  is  N in  each  dimension,  where  D=Nh.  Then 

F(»'j  v')  ~ ^ X 

X-  O ' • ! 

t>*  and  v'  are  each  periodic,  with  period  D.  Values  oft'  and  ^ 

may  be  obtained  at  discrete  values  ofJ  •=  J.  . 

o Nh 

nf  , i) ? Z 'iV  A WX,  oV)  exr>D”r  (**{ 
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where  7°, 


- b Z> 


H-l 


Define  - F(  A , -Z.') 

Azt  = hT) 


In  addition,  for  better  accuracy,  we  introduce  a weighting  factor 
1^x3"  , such  as  that  given  by  Simp son’s  or  the  Trapezoidal  rule. 
Then 


n -1  w-/ 


U/rr  Axcr  + CT 


MT£:o  cr,o 

In  the  same  manner  • <»  / M/ 

f.  _ & p^o‘t-v‘')  .■  In  Equation  (9),  values  of  6< 

represent  positive  values  of  v ' , whereas  values  of  f * ,A/z. 
represent  negative  values  of  [/ 1 . Then 


I2I L.rfef'V/?  * 

(? 

where  A T --  ®, 


* /£? 


0 z/  - • 


NEDMANN-PIERSON  SPECTRUM 


For  a fully  aroused  sea,  the  Neumann-Pierson  Spectrum  is  ^ ^ 

~ TT  & ±- 


— o ■ 

From  References  (l),  (2)  or  (3) 


-rr  5 0 £ rr/t 


*/z. 


whsre  >. 

« =CtCkS*-k+),/i)  .. 

where is  the  argument  of  the  complex  value  k»+,ky.  Then  from 
Equations  (11)  and  (12),  ■ 


5 c** 


- rr  £ 0 < rc 
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= 0,  otherwise 
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FFT  EVALUATION 

Equations  (9),  ( 10)  and  (13)  are  the  main  equations  of  this 
jftudy.  Letting  andky-^  in  Equation  (12)  and  then  defining 

~ Axt  into  Equation  ^9),*  we  may  compute  As 

noted  earlier,  S*  is  the  DFT  of  faj-.  Thus  a two  dimensional  FFT 
may  be  applied. 

One  example  of  the  FFT  will  now  be  discussed.  A wind  speed 
of  5 knots  was  chosen,  h was  set  equal  to^v>-  . This  gives  ^ 
av.  - At-  - _L-  • N was  chosen  as  128. 

Thus  the  range  of  kn  r L -=  S'  — _L 

o 6v/i- 

The  corresponding  /,w  . y ^ ^ Vy  a-  c>  V f jec*' 

and  range  of  ^ ^ <y  ( Xw„t  x <2_&x  8)'^  Y Sec'' 

The  Neumann-Pierscn  spectrum  changes  rapidly  at  the  lower  angular 
frequencies  and  decays  slowly  at  'the  higher  angular  frequencies. 

The  above  choice  of  Akxis.  sufficient  to  samplefine  enough  at  the 
lower  frequencies  &ndMhK„is  sufficiently  large  to  obtain  most 
of  the  contributing  range  of  angular  frequencies. 

The  computer  program  is  shown  in  Appendix  A.  The  first  part 
of  the  program  computes  Equation  (13 ) and  fills  up  an  array  of 
128  x 128  dimensions.  If  each  computation  of  Equation  (13)  is 
considered  as  a complex  quantity,  the  number  of  computer  word 
needed  is  256  x 128  = 32,768.  This  is  within  the  memory  capacity 
of  the  UNIVAC  1108.  If  N=256,  then  the  memory  capacity  must  equal 
256  x 512  = 131,072  words,  to  store  the  complex  values.  This  is 
greater  than  the  capacity  of  the  UNIVAC  1103.  The  elements 
of  the  matrix  were  weighted  byW^j,  where  in  this  case  the  trapezoidal 
rule  was  used.  This  assigns  a on  the  borders  of  the  grid 

and  .25  on  the  corners  of  the  grid  and  J for  all  other  values. 

The  FFT  subroutine  is  then  called.  The  subroutine  is  described 
in  Reference  (4)#.  Following  Equation  (9) , fy,  is  obtained. 

Note  that  d w1  r — =3io-.  and  & (*  n l<)-=  a rr*it- . Thus 

the  distance  between  adjacent  elements  of  the  matrix  = &>-».  . 
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RESULTS 

The  spatial  correlation  at  the  origin  should  be  / . Without 
the  trapezoidal  rule  we  obtain  a value  of  1.04  and  thus  a 4%  error. 
Using  the  trapezoidal  rule,  we  obtain  a value  of  ~ .944  and  thus 

an  error  of  .6$.  The  error  may  be  further  reduced  by  using  the 
collapsed  function  technique.  (see  References  (5)  and  (6)) 

For  a 384  x 384  matrix  collapsed  into  a 128  x 128  matrix,  the 
maximum  value  = .9995  or  a .05$  error.  The  computer  program 
is  written  so  as  to  take  advantage  of  the  collapsing  technique 
-if  desired. 

DeBoer  (Reference  (?))  has  devised  a method  of  obtaining 
the  spatial  correlation  for  the  case  of  the  Neumann-Pierson  spectrum 
as  given  in  Equation  (13)  that  is  more  efficient  than  the  simple 
double  sum.  However,  his  method  may  only  be  used  for  very  special 
cases.  Notiag'that  the  Neumann-Pierson  spectrum  are  a one 
parameter  family,  namely  wind  speed  ^ , DeBoer  multiplies  the 
actual  distance  of  the  spatial  correlation  curve  by1^*.  Following 
-DeBoer,  our  normalized  spatial  correlation  computations  are  presented  • 
-in  Figure  2 -for  the  down  wind  case  and  in  Figure  3 for  the  cross  wind 
case.  -These  curves  agree  with  those  of  DeBoer  and  thus  serve 
is  a check  on  the -FFT  method  and  the  computer  program.  The  correlation 
-for  any  wind  speed  may  be  obtained  by  dividing  the  normalized 
distance  by'1  ^3-  . Note  that  the  correlation  in  the  cross  wind 
direction  is  wider  than  -in  the  dour,  wind  direction.  Note  that  as 
the  wind  speed  -increases-, . the  oc^an  surface  becomes  more  correlated. 

As  the  wind  speed  increases  , the  maximum  peak  of  the  Neumann- 
Pierson  spectrum  shifts  to  lower  frequencies.  These  lower  frequencies 
result  in  a larger  correlation  distance. 

The  Neumanrr-Pi'erson  spectrum  is  symmetric  in  ^ as  well  as 
-in  v1  and  -is  thus  a special  case  of  Equation  (5). 


SUMMARY  AND  CONCLUSIONS 

The  FFT  method  may  be  applied  to  the  directional  wave  spectrum 
to  obtain  the  spatial  correlation  of  the  surface  height.  The 
FFT  method  is  substantially  faster  than  the  straight  forward 
integration  method.  This  work  is  the  first  known  FFT  transform  of 
the  directional  wave  spectrum.  j 


I 


\ r t 


r Tech  Memo 

I No.  T All-120-71 


The  spatial  correlation  values  obtained  by  the  FFT  are 
in  excellent  agreement  with  those  of  DeBoer.  The  curves  presented 
are  dimensionless  - i.e.  they  are  valid  for  all  wind  speeds. 

The  computer  program  and  analysis  in  this  memorandum  should  serve 
as  a guide  to  future  programs  for  the  transformation  of  experimental 
directional  wave  spectra. 
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C SPATIaL  CORKt-LATlOrj  o F SUHFACfc.  

C HEIGHT  VIa  |'HE  KFT 

Phi<Ai«ic.TLi\  NS=x2d»  U = x»*'lDl=l4,i/.ii»NLJlMl=NDl-lrUlMlzNl-l 

PaHA^lT la  ':C— lB  t "IC  1 -!'•>» +-1 , KC2=i'«C+2 

Ii.>JLiCil  JUUOcE  Pfti.CI»l0NlA-H»0-Z)  L 

CUrtPLfcX  a ^ i JL> x r Nj I ) 

DltfcNblO.*  K'N (2)  .r 

OmTh  ivJ/HUl  »HoI/ 

Kc-aL  Al  t JL2)  • 1 1 HOP  ) r uUFFlR ( 1 JOUO) 

Call  PLO  fL,  (bbpFEK » 1 00  JO  ,6) 

Call  PLorio, u.ioO)  

call  plot  to. »o. »-3) 

j1  Fok.'*„1  (bJa.i»uy.4»D5.  Jl 

bo  Fortv.^T  ( 1 a , 214  , 2. jlo,8)  « 

34  Foi<i'iaT  i iAr  4 (2l4,Elo,b)  ) 

CL  — * J 0 ♦ *2 

Cllc-lC-»wL  — — 

□NL=i.D0/LCCC 

Cr.-:=LUo00.J0 

G^JdO.ououO 

GL-0*b  . 

Pi -3#  lAloG2obj56y /9o2*+D0 
o0  Pla£i  (3»3i,lL;.;LJ;:32)  II  , rti  » TS1  r Hbx  , HS2  r 3K , PL 
S^aUxlobL JO .b J/3o JO , JO  . 

S«i-S*:a ^ 


bu=S/ l2.*G) 

SLLm‘jL*bL 

SL  ♦=^>C2«^C2 

oC  4*oC 

PxJ2^Pl/2. 

SiJIuc.=bOKT  £P i j2)  

Va^lF' ‘♦PiuS*  jPIU2*3.J0*SC5 

SibmiUA  T l v AH ) 

UF=CC/WS 

F i »=? , *o/S2 

□o  i 

_DO  2 ub=0/niDiM 

AbJt'i—  0 . 

L>o  3 IP=3fNIMi  

DO  4 vjG=0  ,N1M1 
li<P=li<FlP*Nl)i  __ _ 

IF  (IrvP.Eu.O.AnD* I bC . Eu. 0 ) bGrlS=0 *. 
IF  ( ir<P » Lvn  ( 0 . AijLj  # IS^, , to  • 0 ) GO  To  4 

Il-IPP**2  __ 

Jo=Ibo**2 


RSTSVMABli  COPY 


.__j 


Tech  Memo 
No.  TA11-120-71 


" OiTtG-SGKT  ISwK  j (il/^CLC+JJ/CCCc/^SORnG) 
U,  't.02-0(”C-u**2 

aP  i -2 , *02/  ( 0,  cG2.*S2  ) ) 

Sbnj-^bHb/0MLu**9 

SSi«S=LoHj*G2*LFM/2,00 

Xrs— XKh'/tC 

Yis=I^u/CC 

Tr.£T=«rrti^(YN»XK) 

■ 

AbO?  .^iSiio+^bGH 

4 Coi«T  if,<jE 
3 Co<<Tii'»G)t 

it=ik+i 

j i -ji+ i ■ 

IF  ( 1 K • oL • 0 • OR  • IR^Evv^iv^X'-il ) A jo*'C #oL*G*MSUM 

Ir  ( JS.Eo.  0 .Oi<.  JS»  tO.riul.  ’.I ) ASui'>=  . buO*ASUM 
A 1 1 T , uT ) -ASjJrt 
2.  CuuTinOE 

1 CONTI  

CAwL  FOUa T ( A » ! jU » 2 r — 1 » 0 » 0 ) 


1 1 ( J)  =ANt/ 

iin X Tc.  ( 4 » 3 j ) I * j»y»A|yy ' 

b CONTXivJE 

7 CO.NTIoJt _ 

XI ( NCI ) -0 • 

xnr,C2)=io. _: 

ri((^D=-.6 

ii  u,L2)  = .4 : : 

CALL  L I *^c_  ( X X » Y 1 f NC 1 1 » 0 » 0 ) 

ChlL  AAlolC.  ►O.  »UFNO^MALlZtU  0I3TANCE , -19 1 4 . t 0 . , XllNtl ) r Xl  (NC2)  t 1 
LO.) 

CAuL  AXIt»lO.»0.» HHCC,<RELATl0^f+llr‘+.f90.#Yl(NCl)  »U(NC2)  »10.) 

CALL  Pi-GI  1 20 • # 0 , » -3 ) 

Cau PlO L.lO • i G»  t _953J l 

STOP 


